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Abstract 

Variational wave function is proposed to describe electronic properties of an array of one- 
dimensional conductors coupled by transverse hopping and interaction. For weak or intermediate 
in-chain interaction the wave function has the following structure: Tomonaga-Luttinger bosons 
with momentum higher then some variational quantity A are in their ground state while other 
bosons (with \k\ < A) form kinks - fermion-like excitations of the Tomonaga-Luttinger boson field. 
Nature of the ground state for this quasiparticles can be determined by solving three dimensional 
effective hamiltonian. Since the anisotropy of the effective hamiltonian is small the use of the mean 
field theory is justified. For repulsive interaction possible phases are density wave and p-wave super- 
conductivity. Our method allows us to calculate the low-energy part of different electronic Green's 
functions. In order to do that it is enough to apply standard perturbation theory technique to the 
effective hamiltonian. When the in-chain interaction is strong A vanishes and no fermionic exci- 
tation is present in the system. In this regime the dynamics is described by transversally coupled 
Tomonaga-Luttinger bosons. 
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I. INTRODUCTION 



The adequate description of quasi-one-dimensional (Q1D) conductors remains an unre- 
solved theoretical challenge. Experimentally, at low temperature such systems either three 
dimensional anisotropic Fermi liquids or they freeze into a three dimensional phase with 
broken symmetry[l]. At high temperature their transport properties show many unusual 
features generally attributed to the one-dimensional electron anisotropy. This cross-over 
from ID to 3D is the core problem of Q1D physics. 

It is possible to look at the issue of the dimensional cross-over from another angle. At high 
temperature the proper elementary excitations of the system are Tomonaga-Luttinger (TL) 
bosons. When the temperature is low and the interaction is weak enough the elementary 
excitations are fermions. Therefore, to describe the system at different energy scales one 
needs to explain how high-energy bosons 'cross over' into low-energy fermions. Obviously, 
this is a non-trivial task. 

In this paper we develop a variational approach which accomplishes this goal. To explain 
the structure of the variational wave function let us first consider a one-dimensional con- 
ductor described by TL hamiltonian. The ground state of this system is the ground state 
of TL bosons with all momenta k. Let's turn the transverse hopping on and couple N± of 
these conductors into 3D array. In this situation the system will attempt to lower its ground 
state energy even further by taking advantage of the transverse hopping energy. However, in 
order to participate in hopping the bosons have to form many-body fermion-like excitations 
which have finite overlap with the physical fermion. 

To accommodate for possibility of having two types of excitations, bosonic and fermionic, 
we device our variational state in the following fashion. We introduce intermediate cut-off 
A < A, where A is the cut-off of the ID hamiltonian. All TL bosons whose energy and 
momenta are high (\k\ > A) remain in their ground states. The small momenta bosons 
(| A; | < A) form fermion-like excitations which are delocalized in transverse direction. To 
distinguish between the physical electrons and these fermionic excitations we will refer to 
the latter as quasiparticles. In other words, the wave function can be factorized into two 
parts. High-energy part corresponds to the ground state of \k\ > A TL bosons, low-energy 
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part corresponds to the 3D anisotropic Fermi liquid composed of the quasiparticles. 

The variational energy is minimized by adjusting A. The energy of quasiparticle trans- 
verse hopping is decreasing function of A. At the same time, the in-chain energy grows when 
A grows. The trade-off between the transverse kinetic energy and the in-chain potential en- 
ergy determines the value of A. 

If the optimal value of A is non-zero the low-energy excitations of the system are the 
quasiparticles. Properties of the fermionic quasiparticle state depend on quasiparticle effec- 
tive hamiltonian. It arises naturally after high-energy bosons are 'integrated out'. In this 
effective hamiltonian the anisotropy is insignificant. Standard many-body techniques such 
as perturbation theory and mean field theory can be used to calculate Green's functions 
and map out the quasiparticle phase diagram. Since the physical electron and the quasi- 
particle have finite overlap there is a direct correspondence between the broken symmetry 
phases of the effective hamiltonian and the physical system. We will show that possible 
phases for spinless Q1D electrons with repulsion are the charge density wave (CDW) and 
the superconductivity with the Cooper pairs formed of the electrons on neighboring ID 
chains. 

As the in-chain interaction grows parameter A approaches zero. When A vanish the 
fermionic excitations cease to exist. The system is described by 3D TL boson state. In such 
a regime the ground state is CDW. 

Our approach allows us to obtain some new analytical results. With the help of the 
method it is possible to derive a formula for quasiparticle damping near Fermi surface. 
Also, we evaluate transition temperatures for CDW and superconductivity. The knowledge 
of these temperatures allows us to map out the phase diagram of our system. Although, 
these quantities have been obtained using different numerical techniques [2, 3] the analytical 
expressions had not been reported. 

The paper is organized as follows. In Section II we determine A and derive the effec- 
tive hamiltonian for the fermions. Section III contains the evaluation of the single-particle 
Green's function. Different phases of the effective hamiltonian (and the physical system) 
are mapped in Section IV. The regime where A = is discussed in Section V. We give our 
conclusions in Section VI. 
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II. VARIATIONAL PROCEDURE 



We start our analysis by writing down the hamiltonian for the array of coupled ID 
conductors: 

H = j L dxH, (1) 

w = E«." + E«y> (2) 
nf = iv F (VLvVu - Av^u) + g^u^nv (3) 

n ij = ^ ~ J) £ {^Pi%3 + h ' C -) + ~ j) (^li^Ri^Rj^Lj + h - c -) ( 4 ) 

p=L,R 

+9o(i - j) {$A*l>u + ^Ri^Ri) (VlA; + ^Rj^Rj) ' 

with the real-space cut-off a = ir/A. The fermionic field ip^ creates physical electron with 
the chirality p = L(+) or p = R(— ) on chain i. Transverse interaction constants go (forward 
scattering) and #2% (exchange) are positive. The terms proportional to g and g2k F account 
for the Coulomb repulsion of the electrons on different chains. It is further assumed that: 

9>9o>92k F - (5) 

Now we use Abelian bosonization prescription [4]: 

^l(x) = (2na)- 1 / 2 r ]p e iV ^^ = (2na)' l/2 r lp e^ e ^ + ^\ (6) 

to express the electron hamiltonian in terms of bosonic fields. In the above formula i] p 
are Klein factors, is the TL boson field, <3> is the dual field. The bosonized one-chain 
hamiltonian is: 

H ld [0, $] = ^ (: (V0) 2 : +: (V$) 2 :) + -?-(: (V$) 2 : -: (V0) 2 : ) . (7) 

2 V / 47T v ' 

The symbol : . . . : denotes normal ordering of TL boson operators with respect to non- 
interacting (g = 0) ground state. 

Let us introduce our main variational parameter A < A and use it to split TL boson 
fields into fast (A > — fc F | > A, subscript '>') and slow (\\k\\ \ — kp\ < A, subscript '<') 
modes: 

n ld [©, $] = n ld + n ld = n ld [©<, $<] + n ld [©>, $>] . (8) 
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We define the fermionic field ^ P (x) with the help of equation (6) in which a is substituted 
by a = it/ A and 0< and $< are placed instead of and $. The field \I> is our quasiparticle 
discussed in Introduction. Using this field we re-fermionize H 1 ^. The result is the same as 
(3) with ty p instead of ip p . The transverse terms (4) can be easily re- written if one observe 
that the physical fermion is simply: 

= v /^^t e i v ^(e>+ P *>) ? (9) 

and that the fermionic and bosonic parts in this definition commute with each other. There- 
fore, Hjj is equal to 

Ufj = -(a/a)t(i-j) E (^ P /^« e >^ >^ +p ^-*>^) + h.c.) (10) 

p=L,R 

+(a/a) 2 g 2kF (i - j) (^Ri^M^** - *^ + ^) 

+g (i-j) ((*L*u + *L*Ri) (*L* Lj + + ^V^W^) . 

Our variational wave function has the form: 

|Var) = ni0>;> = |{^» II ^l/vr/C) 1 / 2 exp {-|A;| |$, fc | 2 /K] . (11) 

3 j,k>A 

It is a product of some many-body state ({^pi}) composed of the quasiparticles and the 
ground states | ()>_,■) of ft ld [0>j, $ >i ]. 

Variational ground state energy is found by minimizing the expression: 

£ v = AT ± L^ (A 2 - A 2 ) + ({^ (12) 
W s = £ it* (*L^*u " *L^ W ) + 0*L*u*L*Ri (13) 

i 

ij p=L,R ij 

f = C^, ^ F = C^W C = A/A- (14) 

The number of chains is N±. The TL liquid parameter /C, the electron anomalous dimension 
9 and boson velocity u are defined in the usual way: 

lC = J 2 ^ + 9 > 0= l -(K + K- l -2), u = ±-^2nv F -g) (27rv F + g). (15) 
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The first term of (12) has purely one-dimensional origin. The second term is the energy of 
the quasiparticle ground state. 

Observe that the parameters of the effective hamiltonian t/A and g 2 k F are connected to 
the corresponding bare parameters as if they are subject to the renormalization group (RG) 
flow in the vicinity of the TL fixed point. The explanation to this fact is quite obvious: our 
method of deriving the effective hamiltonian is equivalent to the tree level RG scaling near 
TL fixed point. 

If the transversal interactions are small (g and g 2 k F both less then t/A) they can be 
neglected. In addition, we neglect corrections to the energy due to spontaneous symmetry 
breaking. The latter assumption works when <C 1. Its validity away from this point will be 
discussed at the end of Section V. Under these two conditions the expression (12) becomes: 

eV(ln ± ) « ^ (c 2 - i) a 2 - — e e e mf ■ (i6) 

This variational energy attains its minimum at 

c= f(8P/««pA») 1 ^ if*<l, 
\ if 9 > 1, 

? = E(^)) 2 - as) 

i 

We have to remember, however, that value of the numerical coefficient in (17) is not accurate. 
This is due to the fact that the second term in (16) is calculated under assumption t < A. 
When A gets smaller the coefficient in front of this term acquires some A dependence. We 
neglect the corrections due to this dependence since they are less singular (at small A) then 
the second term of (16). These corrections modify the result for £ quantitatively, therefore, 
it is more appropriate to write 

/ t \V(i-0) 

Coc -) , (19) 



\uA, 

u = y/uvp. (20) 

Using this formula it is easy to show that: 

tocuAocti—j . (21) 

This means that for 9 < 1 the effective transverse hopping amplitude t of the quasiparticle 
^ is of the same order as the quasiparticle longitudinal cut-off energy t> F A. Therefore, due 
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to small anisotropy, the hamiltonian for the quasiparticles (13) can be treated within the 
framework of usual mean field theory and perturbation theory. 

Our calculations, in agreement with renormalization group analysis [5, 6], show that 
for 6 < 1 there is the cross-over energy scale t above which the system is equivalent to a 
collection of decoupled chains while below the transverse hopping becomes important. 

Depending on the interaction and the anisotropy the region 6 < 1 can be further split 
into two parts. The transverse hopping contribution to the variational energy (second term 
of (16)) can be re- written as follows: 



If the argument of the exponential function is small the exponential can be replaced by the 
first few terms of the Taylor series. In such a situation the contribution of the in-chain 
interaction to the total energy, eq. (16), can be calculated perturbatively. One-dimensional 
effects are virtually unobservable. This is the weak coupling regime. 

When the anisotropy and the in-chain interaction are strong the exponential cannot be 
approximated accurately by the low order Taylor expansion. The system is in the inter- 
mediate coupling regime now. In order to obtain a reliable answer in such a regime it is 
not enough to apply finite-order perturbation theory. Our method converts the system of 
physical electrons with intermediate coupling into the system of quasiparticles with weak 
coupling. The latter can be studied by standard perturbation theory. 

As a function of the bare transverse hopping amplitude t the cross-over from the weak 
coupling to the intermediate coupling occurs at: 



For the weak in-chain interaction <C 1. If this is the case it is necessary to have exponen- 
tially small transverse hopping amplitude t in order to observe non-trivial Q1D effects. 

When > 1 the effective cut-off momentum A is zero. The quasiparticles are not formed. 
The system can be viewed as a collection of TL bosons weakly coupled by the transverse 
exchange interaction. The possibility of such state was first pointed out by Wen [7]. It is 
natural to call such a regime strong coupling. Section V is reserved for discussion of strong 
coupling. 




(22) 




(23) 
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III. SINGLE-ELECTRON GREEN'S FUNCTION 



The calculation of different propagators for Q1D system is an open question. Our ap- 
proach allows for easy evaluation of the low-energy part of Green's functions in the inter- 
mediate coupling regime. The high-energy parts of Q1D Green's functions are believed to 
coincide with the Green's functions of TL model. The latter have been discussed extensively 
in the literature. 

The Matsubara propagator of the physical electronic field ifa is equal to: 

G L (x, Rj., r) = I (T {* Li (x, r) ^ (0, 0)})^ x (24) 

IT | e - i v / ^(e> l ( :E ,T)+$ > ,( :E ,T)) e i v ^(0 >J (o,o)+<i> >J (o,o))|\ ) 

R._l = Rj — Rj. (25) 

The notation (. . .)^ stands for averaging with respect to the quasiparticle ground state |0$). 
Likewise, (. . ) > stands for the expectation value with respect to |0>) state. 
The bosonic part of this formula can be immediately calculated: 

1 (j | e -iV5F(e> i ( a!l r)+*> i (x ) r)) e iVi(e> i (o 1 o)+*> i (o 1 o)) = ^Gj^j s .. + _ 5 ,,y (26) 

Here Ql d (Gl d ) is the Matsubara Green's function of the Tomonaga-Luttinger model with 
the cut-off A (A). 

Our variational wave function does not take into account correlations between $>j(0>j) 
and $>j(0>j) if i ^ j. However, the above formula is correct at least for large r > 1/uA 
or small frequency u < uA where those correlations are not important. In such a limit the 
boson part of (24) is a constant equal to ( 9 . 

Once the bosonic propagator is found it is necessary to calculate the quasiparticle Green's 
function. This can be done with the help of standard diagrammatic technique. If we neglect 
interactions between the quasiparticles the single-electron Green's function is equal to: 

G L (iu;,p||,p) = — —2 -, (27) 

V > iuj + v F p\\ - 

where the renormalized transverse kinetic energy is given by: 

4 = - 2 C e E^) C0S (P- R *)- ( 28 ) 

i 

This result coincides with the Green's function derived by RG [2, 6]. 



Our method allows to improve the above formula for the single-electron propagator by tak- 
ing interaction between the quasiparticles into account. Neglecting (i) symmetry-breaking 
which becomes important for very small frequency only and (ii) the transverse couplings g 
and c/2k F (see (5)) one can identify three second-order diagrams contributing to the single- 
quasiparticle self-energy (fig. 1). They are: (a) scattering on the polarization bubble of the 
same chirality as the incoming quasiparticle, (b) scattering on the polarization bubble of the 
opposite chirality and (c) the vertex correction. The diagrams (a) and (c) are identical in 
magnitude and opposite in sign. Thus, (b) is the only diagram on fig.l which needs to be 
evaluated. 

First, we calculate the quasiparticle polarization bubble Vr. It equals to: 

v» (in, h , k) = / aiq 5 (? m + $ (29) 

where the notation / . . . = (27r)~ 3 5 2 / dq\\d 2 q. . . is used. The symbol b denotes the trans- 
verse lattice constant. The self-energy is equal to: 

S L (iu;,p||,p) =-T£ / g L (iuj + iQ, Pll +k h p + k)v R (in,k h k). (30) 
n Jk \\ k 

After summing over Q the following expression for the self-energy is derived: 

E L = -«? 2 / / (^,1+^-^)5(^11+5^) x (31) 

«/ fen k Jgyq 

(v F k\\ + £^ +q + 4) - n F (v F (k\\ + p||) + e£ +p ) 
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iu + 2v F k ]l + vfP\\ - £q- + £"£ +q - e£ +p 
(v F k\\ + £^ +q - £q-) 5 (v F q\\ + £- 



+ 



g I > 



fc||k Ayq sinh { (w F /c|| + e^ +q + e^j /t} (iu; + 2v F k\\ + ^ F P|| - £ q + ^k+q - e k+ Py l 

When T = the first term can be further simplified. The Fermi distribution n F becomes 
the step-function. In such a situation it is possible to perform integration over k\\ and p\\ 
exactly. The second integral in the above equation appears due to the relation between 
the Fermi distribution n F and the Bose distribution n^,: n^(uj) + n F (uj) = 1/ sinh(a;/T). 
At zero temperature this integral vanishes. In the resultant T = expression for S L the 
transverse kinetic energy e 1 - always enter in the combination £ q+k + £k+ p ~~ £q- Therefore, 
it is convenient to introduce the quantity: 

^(e\ P) = / d ^bH {^ - eU - 4 +p + si) . (32) 



With this definition the self-energy can be compactly written as follows: 

S L (iu,p h p) = — p|| (33) 

9 2 f , j_ i / i s /. l 4A 2 

-T6^f7 d£ v {e 'P)(--^n- £ ) lo ^ 2 + ( ^ ||+£ , )2 - 

The Green's function of the physical electron is equal to ( 6 (iu; + v-pp\\ — — S L ) _1 . Note, 
that the logarithmic divergence of the self-energy a hallmark of the Fermi liquid picture 
break-down in TL model, is capped in the presence of the transverse hopping. This justifies 
the use of the perturbation theory. 

By doing analytical continuation of (33) it is possible to calculate the retarded self-energy 
E ret whose imaginary part is the quasiparticle damping: 

7 = -ImS ret = £^^(-v m ,p)u; 2 . (34) 

The transverse density of states can be estimated as u 1 - oc 1/t. This gives us 7 oc (g /vf) 2 uj 2 ft. 
On the mass shell uj = — v-pp\\ + £ p the expression for 7 becomes: 




where (py — p-p) is the distance from a given point (py, p) of the Brillouin zone to the Fermi 
surface ^fPf = £p along x direction. 

We need to issue a warning in connection to the accuracy of El- It is not correct to think of 
(33) as 0(g 2 ) expression for the physical electron self-energy. Indeed, the physical electron 
Green's function (27) already contains all orders of g entering though the quasiparticle 
renormalization ( and renormalized transverse hopping e. It is necessary to remember 
that our variational approach is uncontrollable approximation. It lacks a small parameter 
controlling the quality of the results. Therefore, it is not clear how accurate the expression 
(33) is. 

In Ref. [2] the self-energy was evaluated numerically for the system with infinite trans- 
verse dimensions. However, those calculations are more complicated technically and do not 
give analytical answer for the self-energy. 

IV. PHASE DIAGRAM 

In Section II we derived the low-energy effective hamiltonian for the quasiparticles. Now 
we will apply the mean field theory to obtain the phase diagram of the effective hamiltonian. 
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The experimentally observable phase diagram for the physical electrons coincides exactly 
with that of the quasiparticles. To prove this let us calculate (V'uV'Ri/ f° r T ^ v-A: 



The physical CDW order parameter is proportional to the CDW expectation value of the 
quasiparticles. Similar formulas can be obtained for other order parameters. For example, 

(ViA) = C 1 ^" 1 (*Li*L)<p and (WA) = C 6 (*L*!y)*- Therefore, we can determine 
the phase diagram of (1) by mapping the phases of the hamiltonian (13). 
We consider four order parameters. One is the charge-density wave: 

p 2kFl = (37) 

and there are three types of the superconducting order: 

A±y = \ (*l*h ± , (38) 

Ac, = *L*L- (39) 
The in-chain potential energy can be re-written in terms of p and A in the following manner: 

^*L*Li*Ri*Ri = -9P2k F Ak F i = g&oAoi- (40) 

The exchange interaction can be expressed as: 

~92k F (* L*Ri*iy*Li + h.c.) = ~g 2 k F {p2k F Ak Fj + h -c-) = 2g 2kp (A_, j A t _ ij - A +y .A^) . 

(41) 

Finally, a part of the transverse forward scattering which describes the interaction between 
the fermions of different chiralities is equal to: 

9o (*L*iA*r, + *L*rA-* Lj -) = 2^0 (A +iJ A^. + A_yAL y ) . (42) 

The part of the forward scattering which accounts for the interaction between the fermions 
of the same chirality cannot be expressed in terms of these four order parameters. 

The effective coupling for CDW is always bigger then the effective coupling for the su- 
perconducting order parameter A + : 

gcuw > g sc , where (43) 
fl-CDW = 9 + z±hk F , (44) 
#sc = 92k F ~ 9o = C 21C ~ 2 92k F - go, (45) 
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and z± is the coordination number for a chain. Thus, at T = for the perfect nesting the 
system is always in CDW phase with A + order parameter phase being meta-stable (g sc > 0) 
or unstable (g sc < 0). Other order parameters, A and A_, are unstable. 

When the external pressure is applied the amplitude t 2 for hopping to the next-to-nearest 
chain begins to grow and spoils the Fermi surface nesting. This undermines stability of CDW 
and drives the transition temperature to zero [8] . Indeed, in the latter reference the following 
simple estimate for the density wave susceptibility was obtained: 



1 log (2v F A/T , if T>t 2 = (% , 

X oc x < ) _ ; _ (46) 

2irv F \ log (2v F A/t 2 ) , HT<t 2 = (%, 

The CDW transition temperature is derived by equating (g + z±g 2kp )x an d unity. For t 2 = 
it is: 

t cdw K w fA exp (-2nv F / (g + z ± g 2kF )) . (47) 

If t 2 > the transition temperature Tcdw becomes smaller then T c ^ w . It vanishes when 
t 2 oc T C d W . That is, exponentially small t 2 is enough to destroy CDW. 

What happens after CDW is destroyed depends on the sign of g sc . If g sc > the 
ground state is superconducting. Otherwise, it is the Fermi liquid. We can perform the 
same type of analysis we did above for CDW. The superconductivity is rather insensi- 
tive to the nesting properties of the Fermi surface. The susceptibility for A + is equal to 
(l/27rcTO F )log (2v F A/T), where a is a constant of order of unity. The critical temperature 
is found to be: 

T c oc v F A exp (-2nav F /g sc ) , (48) 

if g sc > 0. Even when g 2 fc F < 9o the effective coupling g sc may be positive provided that 
the in-chain interaction is repulsive (JC < 1) and the electron hopping anisotropy parameter 
(uA/t) is big: 

On /? -yA\ (2-2K)/(l-0) 



C^- 2 >^^ _ > > I. (49) 

92k F V t J g 2kp 

For the system in the intermediate coupling regime this condition is likely to be satisfied. 

It is interesting to note that the external pressure detriments not only CDW but the 
superconductivity as well. Under pressure the anisotropy parameter (uA/t) decays. The 
superconducting transition temperature gets smaller as the anisotropy decreases. At pressure 
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higher then some critical value the condition (49) is no longer satisfied. In this region the 
superconductivity is unstable and the ground state is the Fermi liquid. 

The qualitative phase diagram is presented on the fig. 2. It shares two remarkable fea- 
tures with the phase diagram of the organic Q1D superconductors [1]: (i) the density wave 
phase and superconductivity have common boundary; (ii) the superconducting transition 
temperature vanishes at high pressure. 

Our order parameter A + deviates from the more common version A . The order pa- 
rameter A + was proposed quite some time ago [9]. Recently, this suggestion found further 
support in the renormalization group calculations of Ref. [3]. The advantage of A + stems 
from the fact that by having two electrons of a Cooper pair on different chains we avoid 
increasing in-chain potential energy. 

The origin of the superconducting phase in our system is an interesting question worth 
discussing in more details. In conventional BCS model the superconductivity is stable be- 
cause it minimizes the potential energy of the electron-electron interaction. We can make 
this claim rigorous by considering the following derivation. BCS hamiltonian density 

2m 



n BCS = r + v = 5>l(!--*'W,- #}vl^ T (50) 



consists of two term: kinetic energy density T and potential energy density V. At zero 
temperature the superconducting state energy density £ sc = (H) sc is smaller then the normal 
energy density £ n = (Ti) n . This condensation energy density 

£ c = £ Q - £ sc oc z/T c 2 , (51) 
v = ix~ 2 m 2 v-p, vf = ^2/i/m, (52) 

is entirely due to depletion of interaction in the superconducting state: 

(V> n - (V) sc > 0. (53) 

As for the kinetic energy it grows in the superconducting state: 

COn " <Uc < 0. (54) 

To prove this we will use Feynman formula which allows to calculate the ground state 
expectation value of any term cO of the hamiltonian density: 

c(0) = c^, (55) 
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S c oc ul^mWexpi-afi-Wm-Wg- 1 ) 



(56) 



(57) 
(58) 



where wd is Debye frequency and a > is a constant of order unity. The inequalities (53) 
and (54) immediately follow from the expressions above. These inequalities mean that it is 
the electron-electron attraction which triggers BCS superconductivity. This fact is a very 
well known fact of the superconductivity mean-field theory. 

However, in the system with strong repulsion, such as Q1D or high-T c materials, it 
is difficult to construct a mean-field superconducting phase which lowers the interaction 
energy. Our model for which we develop the consistent many-body approach can be used to 
discuss this issue beyond the mean-field approximation. 

For our model it is easy to determine that the transverse forward scattering energy is 
increased and the exchange energy is decreased by the superconductivity. This result is a 
direct consequence of (42) and (41). 

Contributions of other terms can be found with the help of Feynman formula. The 
condensation energy density is of the order of — T c 2 /t> F . Thus, differentiating the critical 
temperature (48) with respect to some coupling constant of (1) we can determine how a 
ground state energy contribution of a given term is modified by presence of the supercon- 
ductivity. A derivative of the critical temperature with respect to a parameter x is equal 
to: 



provided that g sc <C v-p. Combining this result with (45) we conclude that in the supercon- 
ducting state the transverse hopping energy is higher: 



(#tA, 4^) n - (*&C*u 4^) sc - 9^ log (±) | (^) > 0, (61) 




(59) 




(60) 



and the in-chain potential energy is lower then in the normal state: 
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since both log(t/wA) and the derivative with respect to g are negative. 

We have proven that in our case the superconductivity is triggered by the electron- 
electron repulsion. This result is quite unexpected. It is has a many-body nature and 
cannot be obtained within a mean- field theory for the hamiltonian (1). This mechanism of 
superconductivity is very similar to the Kohn-Luttinger proposal. Classical Kohn-Luttinger 
mechanism predicts extremely low critical temperature. In our case, however, the effective 
coupling constant g sc is a non-analytical function of the bare parameters. As a consequence, 
our transition temperature (48) does not have to be small. 

V. STRONG COUPLING REGIME 

We have seen above that if > 1 then A is zero. This means that quasiparticles are 
not formed and it is more convenient to treat the system in terms of TL boson only. The 
bosonized hamiltonian (1) has the form: 

n] A = |(/c (ve,) 2 + /c- 1 (v$,) 2 ) , (62) 

In this formula both the transverse hopping term which is irrelevant in RG sense and the 
forward scattering term which is marginal are omitted. Their effect is small as compared 
with that of the strongly relevant exchange interaction, eq. (63). 

The relevance of the exchange interaction indicates that at low temperature the system 
freezes into a state with the finite expectation value 7^ 0. This phase is CDW. It can 
be easily proved by bosonizing CDW order parameter: V'uV'Ri ^ (27ra) _1 exp(iv / 47r$j). The 
finite expectation value of the field $ is inherited by CDW order parameter. 

We describe this regime with the help of our variational wave function. Since A = one 
can write the wave function in terms of TL bosonic field only: 

l Var > = II -j^=exp{-\^ k \ 2 /Aal}. (64) 

k>0,i y 27T<7 fc 

This expression is a slight generalization of (11): in the latter equation the parameters a\ = 
K/A\k\. Here we do not fix a\. Instead, they will be determined variationally. Variational 
energy is: 

E "i {LN ^tfAk^) -^ a2 -pH s< *4 (65) 

15 



Minimizing this energy with respect to a\ we find: 



o\ = . Ku (66) 

4^ 2 fc 2 + a c ; 



*CDW 



a cdw = 8ng 2 k F lCuA 2 exp < -2/C / , £ = > oc g 2 k F uA 2 f A °° W > ) • (67) 

I 7o V £2 + A cdwJ V mA 7 

The quantity Acdw has the meaning of the excitation gap due to CDW order. This gap, 
together with the transition temperature, can be found by solving the last equation: 

(92k V/P-aK) 

T CDW oc A CDW oc uA [—^j ■ (68) 

The variational energy is: 

E W /(LN ± ) (x A 2 CDW /u. (69) 

These results are correct when Ti 1 - couples only those chains which are nearest neighbors. 
The next-to-nearest neighbor coupling frustrates CDW phase. We will not discuss the effect 
of the frustration in this paper. 

Finally, let us discuss cross-over from strong to intermediate coupling regime. Such cross- 
over occurs when the intermediate coupling Fermi liquid energy, eq. (16), becomes equal to 
the strong coupling CDW energy, eq.(69): 

/ f n, 1/(1-0) /„ , x l/(2-2K) 

P/v F oc A CDW / M or ±) oc . (70) 



,uAJ V u 

This equation defines 9 c (t,g 2 k F ) < 1 at which the cross-over takes place. At small 9 < 9 C the 
system behaves as the Fermi liquid whose properties we discussed in the previous sections. 
When 9 > 9 C the expression (17) is no longer applicable: the necessary requirement for 
smallness of the energy associated with symmetry breaking is violated. The expression (69) 
has to be used instead. 

Fig. 3 shows how the strong coupling regime at big 9 is replaced by the intermediate 
coupling regime at smaller 9. The transition temperature of CDW, eq.(68), drops sharply 
and becomes exponentially small, eq.(47), as 9 gets smaller then 9 C . This diagram was 
discussed in [10] for a similar model. 



VI. CONCLUSIONS 



We propose in this paper the variational wave function for Q1D system. Our procedure 
key ingredient is the splitting of TL bosons into high-momentum and low-momentum modes. 
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While high-momentum modes are in their ground state the low-momentum modes form 
quasiparticles which delocalize in the transverse directions. 

Our method can be viewed as a variational implementation of the lowest order RG scaling 
near TL liquid fixed point. When the transverse hopping amplitude becomes of the order 
of uA the scaling must be stopped. The renormalized hamiltonian should be treated as the 
hamiltonian for the quasiparticles. 

Our method gives us a possibility of computing different Green's functions beyond RG 
using standard diagrammatic technique. As an example we calculated the lowest order 
self-energy for the one-particle propagator. 

Depending on the strength of the in-chain interaction and the anisotropy the system 
may be in one of three regimes: strong, intermediate or weak coupling. In the strong 
coupling regime quasiparticles are not formed and the system is better described in terms 
of TL bosons. In weak and intermediate coupling regime the low-lying degrees of freedom 
are quasiparticles. The ground state of these fermions may be either Fermi liquid, the 
superconductivity or CDW. The phase diagram of our Q1D model looks very similar to that 
of the organic Q1D superconductors. 

Unlike classical BCS superconducting phase, the one in our model is stabilized without 
any attraction between the electrons. It is similar to Kohn-Luttinger superconductivity. 
However, our effective coupling constant is bigger then that of Kohn-Luttinger. This guar- 
antees that the critical temperature in our model is not unacceptably small. 
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FIG. 1: Lowest order contribution to the self-energy of the quasiparticle. 
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FIG. 2: Qualitative phase diagram of our model. Solid lines show second-order phase transitions 
into CDW and the superconducting phase. Dashed line shows the first-order transtion between 
CDW and the superconductivity. 
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FIG. 3: The energy scale assosiated with transverse hopping A decreases when 9 grows. The CDW 
transition temperature Tcdw increases as 6 grows. At 9 C where both energy scales are of the same 
order the cross-over from the intermediate to strong coupling occurs. 
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